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^ ! Abstract 
l> 

\^ , We investigate the physical characteristics of the Broad Line Region (BLR) of NGC 5548 using the 

Boltzmann plot (BP) method, that we applied on the Balmer lines observed from 1996 to 2004. We find 
that variability earlier detected in the lines, is also present in the temperature parameter A, obtained from 
IQ ■ the BP method. Using the BP temperature parameter A to calculate temperature, we find that the average 

temperature for the considered period was Tw 10000 K, and that varies from 5000 K (in 2002) to 15000 
K (in 1998). This variation correlates with the AGN-componcnt of the optical continuum flux (r = 0.85) 
Qh' and may indicate the existence of an accretion disc in the BLR of NGC 5548. 

Q . Key words: galaxies: active - galaxies: individual: NGC 5548 - galaxies: Seyfert 

' 1. Introduction not explain the measured line flux ratios. The classical 

^ studies point toward photoionization, as the main heat- 

The physics and kinematics in the Broad Line Region ing source of the BLR emitting gas, that may explain 

^ , (BLR) of Active Galactic Nuclei (AGN) is more compli- observed spectra of AGN (see e.g. Kwan & Krolik 1981, 

■ cated than in its Narrow Line Region (NLR) or in gaseous Osterbrock 1989, Baldwin et al. 1995, Marziani et al. 

■ ■ ' nebulae (Osterbrock 1989; Krolik 1999; Sulentic et al. 1996, Baldwin et al. 1996, Ferland et al. 1998, Krolik 

2000 and references therein). In contrast to the NLR 1999, Korista & Goad 2004). On the other hand, some 

where forbidden lines (e.g. [OIII] and [Nil] lines) can be authors e.g. Dumont et al. (1998) have given favor to a 

used as emitting plasma diagnostics, the physical condi- non radiatively heated region that contributes to the BLR 

tions in the BLR cannot be obtained using simple relations line spectrum, which may be the case in the BLR of NGC 

between the line ratios. The pure recombination condi- 5548. Therefore, the photoionization, recombination and 

tions cannot be applied in BLRs, e.g. the flux line ratios collisions could be considered as relevant processes that 

arc different than expected in the case of recombination occur in BLRs. At larger ionization parameters, reconibi- 

(e.g. in some AGN Lya/HP w 10, Osterbrock 1989). nation is more important, but at higher temperatures or in 

Several effects can result in such line flux ratios of hy- the case of low ionization parameters the collisional exci- 

drogen lines, particularly the collisional excitation and ex- tation becomes important also (Osterbrock 1989). These 

tinction effects. For over more than 30 years now, nu- two effects, as well as radiative-transfer effects in Balmer 

merous papers have been dedicated to this problem, see lines, should be taken into account when explaining the 

e.g. Netzer (1975), Ferland & Netzer (1979), Ferland et ratios of Hydrogen lines. Moreover, the geometry and pos- 

al. 1979, Kwan (1984), Rces et al. (1989), Ferland et al. sible stratification in the BLR may affect the continuum 

(1992), Shields & Ferland (1993), Dumont et al. (1998), and line spectra of AGN (Goad et al. 1993). 

etc. Dust is present in the host galaxy of an AGN (sec Recently, Popovic (2003, 2006ab) showed that in the 

e.g. Crenshaw et al. 2002, Crenshaw et al. 2004, Gabel BLR of some AGN, the Balmer line ratios follow the 

et al. 2005, etc.), but it seems that in some cases it can- Boltzmann plot (BP), which indicates that the popula- 
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tion of the levels with n > 3 follow the Saha-Boltzmann 
equation. If the population of excited levels of the Balmer 
series can be described by Saha-Boltzmann equation, one 
can determine the temperature of the region where these 
lines are formed. This method, and also recently given 
method by Laor (2006) allow us to ascertain information 
about the physical (thermodynamical) properties of the 
BLR in a direct way by only measuring the broad line 
parameters (Ilic et al. 2007). 

In order to test BP method as a tool for the temperature 
diagnostic in the BLR, we investigate here the spectra of 
NGC 5548 observed from 1996 to 2004 (Shapovalova et al. 
2004). First, we found that BP method can be applied to 
the broad Balmer lines of NGC 5548, and after that we 
determined the temperature parameter A using the BP 
method. 

NGC 5548 is one of the most intensively monitored Sy 1 
galaxies (Shapovalova et al. 2004 and references therein) 
and the physics of its BLR have been studied by many 
authors (e.g. Krohk et al. 1991, Shields & Ferland 1993, 
O'Brien et al. 1994, Dumont et al. 1998, Vestergaard & 
Peterson 2005, Korista & Goad 2000, 2004). In particular, 
Dumont et al. (1998) discussed the physical model of the 
NGC 5548 BLR, and found inconsistencies between the 
photoionization model and the measured line fluxes, line 
flux variations and the underlying continuum. The aim of 
this work is to investigate the variations of the physical 
properties of the NGC 5548 BLR using the BP method 
and Balmer lines, during a period of 8 years. 

2. The BP method 

For plasma of the length £ that emits along the line of 
sight, the flux (or the spectrally integrated emission- line 
intensity) F^i of the transition from upper to lower level 
{u^ I) can be calculated as (see also Popovic 2006b): 



he r 
Fui = —guAui J Nu{x)dx 



(1) 



where A is transition wavelength, g„ statistical weight 
of the upper level, A-^i transition probability, 7V„ is the 
number density of emitters excited in upper level which 
effectively contribute to the line flux (which are not ab- 
sorbed) and h and c are the well known constants (Planck 
and speed of light). In general, the can be inhomo- 
geneous across the line of sight and also, the radiative 
self- absorption can be present. But, assuming that pop- 
ulation in the observed region (in all layers) follows the 
Saha-Botzmann distribution one can write 



Nuix) « ^M^exp(-K/fcr(x)), 



(2) 



where Z is the partition function, Nq is the total number 
density of radiating species, Eu is the energy of the upper 
level, T is temperature and k is the Boltzmann constant. 
In the case of optically thin plasma^ with relatively 



Note here that we cannot expect that plasma in the BLR is opti- 
cal thin, further in the text we will consider possible distribution 
of emitters within the BLR, see Eq. (5) 



small changes in the density and temperature one can 
write (see, e.g., Gricm 1997; Konjevic 1999) 

he A A -^0 

Fui^^guAui J Nudxw—Auigui-^cxp[-Eu/kT).[3) 

For one line series (as e.g. Balmer line series) if the pop- 
ulation of the upper energy states (n > 3)^ adheres to 
a Saha-Boltzmann distribution than one can determine 
their temperature (T), from a Boltzmann plot as 



log(i^„) = log = B- AEu, 



(4) 



' OuAui 

where Fn is normalized flux, B and A arc BP parame- 
ters. Parameter A, deflncd && A ~ logjo(e)/fcr (where 
e « 2.718), is the temperature indicator out of which we 
can estimate the temperature of the region where these 
lines are formed. 

On the other hand, one cannot expect the homoge- 
nous distribution of the physical parameters and density 
of emitters across the line of sight in an extensive BLR. 
But, if we still have that population follows the Saha- 
Boltzmann equation, Eq. (1) can be written as: 



F, 



he 



ul 



-T-9uAul 

A JQ 



Noix) 



cxp{~ Ey^ / kT (x))dx 



(5) 



and if we divide the BLR in small layers with the same 
physical conditions and emitter density, then we can write: 



he 



Fui — —guAui 



A^o(^) 
Z 



eM-Eu/kT{i))£, 



(6) 



If we assume that the temperatures across the BLR vary 
as T{i) = Tav ± ATi and emitter density as No{i) = TVo™ ± 
AAio(«), the Eq. (6) can be written as: 



he N"'^ " ^ 



E,, 



kTa.{l + 5Ti, 



where 5Ti = l^Ti/Tav and 5Na = ANq/NS". If in a BLR 
the temperature and emitter density does not vary too 
much, i.e. ANq/Nq «1 and ATi /Tav << 1, then the 
Eq. (7) becomes 



he N"''" 

-j-9uAui—^eyip{-Eu/kTav)i. 



(8) 



meaning that the Eq. (4) can be used to determine Tav 
in the BLR (see Popovic 2006b). 

Radiative transfer issues (e.g. high optical depths of 
the Hydrogen lines) complicates the picture, but some 
additional tests are possible (see §5.2 and also Ilic 2007). 
Furthermore, the "Case B" recombination line ratios can 
produce similar Boltzmann plots of the Balmer line series 
in some AGN whose BLR is close to " Case B" recombina- 
tion (see Osterbrock 1989). Due to the physical conditions 
(densities and temperatures) in such BLRs the constant 
A is too small {A < 0.2) and the Boltzmann plot cannot 
be applied for diagnostics of the temperature even if the 



^ Note here that since the emission deexcitation goes as m — > Z it 
is not necessary that level I has a Saha-Boltzmann distribution. 



No. ] 

BP can be properly applied (Popovic 2003, 2006ab, Ilic at 
al. 2006). 

We should emphasize that the BP method (Popovic 
2003, 2006ab) docs not take into account any a priori 
physics in the BLR, besides that the Balmer lines are orig- 
inating in the same emitting region. The method includes 
the intensities and transition parameters of five Balmer 
lines, not only the ratio of two or three lines, as usually 
considered. We would like to point out that the BP is not 
the Balmer decrement, i.e. the ratio of the Balmer lines. 
Also, with this method, we concentrate on the thermody- 
namical state of the BLR trying to sec if the BP method 
works at all, without assumptions on any "macroscopic" 
model caused by the "microscopic" physics. 

3. Observations and selection of the spectra 

Optical spectra of NGC 5548 were taken with the 
6 m and 1 m telescopes of SAO (Russia, 1996-2004) 
and at INAOE's 2.1 m telescope at the Guillermo Haro 
Observatory (GHO) at Cananca, Sonora, Mexico (1998- 
2004). They were obtained with a long slit spectrograph 
equipped with CCDs. The typical wavelength range cov- 
ered was from 4000 A to 7500 A, the spectral resolution 
was 4.5-15 A, and the S/N ratio was >50 in the contin- 
uum near the Ha and H/3. Spectrophotometric standard 
stars were observed every night. Details concerning the 
spectroscopic observations are given in Shapovalova et al. 
(2004). 

Note that the spectra of NGC 5548 were obtained un- 
der the AGN spectral monitoring program with the aim 
to study the variability of the broad Balmer emission lines 
Ha and H/3. The errors of the continuum at 5100A and 
line fluxes were about 3-5%. From the spectral database 
collected during the monitoring program of NGC 5548 
(wl50 spectra) we selected only 24 spectra (Table 1), 
which satisfy the following: 

1. good photometric conditions; 

2. S/N > 10 in the continuum near the He line; 

3. spectrum covers the wavelength range from A4000A 
to A7000A; 

4. the broad component of the Balmer lines from Ha 
to He are presented. 

The basic data of these observations are presented in 
Table 1. 

The spectrophotometric data reduction was carried out 
either with the software developed at the SAO RAS by 
Vlasyuk (1993) and with IRAF for the spectra obtained 
in Mexico. The image reduction process included bias 
subtraction, flat-field corrections, cosmic ray removal, 2D 
wavelength linearization, sky spectrum subtraction, stack- 
ing of the spectra for every set-up, and flux calibration 
based on standard star observations. We also remove the 
continuum of the host galaxy. 

The procedure of absolute calibration of the spectra is 
described in details in Shapovalova et al. (2004) and will 
not be repeated here. 
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Fig. 1. The adopted narrow emission line template. 

4. Measurements of the broad line flux 

One can expect that the BP can be applied only in the 
case of the broad Balmer lines (Popovic 2003), and first 
step was to subtract the continuum and the narrow and 
satellite lines from the Balmer lines. To estimate the con- 
tribution of these lines we used a template of narrow and 
satellite lines (Figure 1) estimated from two independent 
measurements in the case of the minimum intensity of the 
broad Balmer lines. In this template the Fe II lines in the 
H/3 spectral region arc included. The optical Fc II lines 
showed variability with smaller amplitude than that of 
the H/3 (around 50%-75% of H/3, Vestegaard & Peterson 
2005). Testing the contribution of the Fe II residuals (af- 
ter subtraction of the template), we concluded that it may 
contribute about 2-5% to the total measured line fiux of 
the H/3 line. 

Furthermore, wc a priori take into account that the con- 
tinuum subtraction contributes to the error-bars within 
10% of the measured fluxes. In general, we can expect 
this value in weak Balmer lines (e.g. He), but in the case 
of the Ha and H/3 lines, it should be even smaller. 

The fluxes of the Balmer lines were measured several 
times, but it is important that independent measurements 
were performed by two persons. Then, the error-bars were 
calculated as: 

AF, = Ai^"'"'' + 0.1-i^r'=" (9) 
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Fig. 2. The Boltzmann Plots of the Balmer line series of NGC 5548 from different periods. _F„ is calculated using the Eq. (4) after 
the normalization of all Balmer lines to the Hji flux (F^i{HI3) = 1). 



where F™<^« is the measured flux, AF™*^* is statistical er- 
ror (within Icr) obtained from several measurements and 
Q ipmes taken to be the error of the continuum and 
narrow line template subtraction. The flux ratios of the 
Balmer lines and the fluxes of the H/3 from different pe- 
riods are given in Table 2. The error-bars presented in 
the Table 2 were obtained as (see e.g. Agekyan 1972, 
Bevington & Robinson 2003): 



Ai?i — Ri 



prnes 



(10) 



where Ri is the ratio of Fi/Fnp- 

In Figure 2 the error-bars {Alogio{Fn)) have been cal- 
culated as: 

AF 

Alogio(F^) = ^ . (11) 

yi>A n, ^„,es X logglO ^ ' 

One can conclude from Eq. (4) that the absolute scale 
is not important for the BP analysis, which uses the slope 
to determine the temperature coefficients (A), but one 
should re-scale fluxes of Balmer lines for the same factor 
as well as the AFi. The parameters A are obtained from 
the best fit of the measured values and we give asymptotic 
standard error AA (within Icr) in Table 2. 

We should note here that we used only those observed 
spectra where the spectral region from He to Ha was cov- 
ered, except for the two spectra (see Table 2) that covered 
only the interval from H^ to Ha. We used the aforemen- 
tioned spectra in order to have the data in years 2000 and 
2001. 

Reddening 

The reddening effect can influence the Balmer line ratio 
(e.g. Crenshaw & Kraemer 2001; Crenshaw et al. 2001; 
2002; Popovic 2003) and consequently the temperature 
parameter obtained by the BP. In the case of NGC 5548 
the Galactic reddening is negligible E(B-V) = 0.020 mag 
(Burstein & Heiles 1982; Schlegel et al. 1998), and here it 
is not considered. Concerning the intrinsic reddening, in 
first approximation we can adopt that it might be 30% - 
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Table 1. The basic data of the spectroscopic observations of NGC 5548. The column "Code" refers to the telescope used for the 
observations: L - 6m telescope SAO RAS (Russia), GH - 2.1m telescope GHO (Cananea, Mexico). 



File 


UT-Date 


JD 


Code 


Aperture 


Sp. range 


Resolution 


Seeing 


(spectra) 




(2400000+) 




(arcsec) 


(A) 


(A) 


(arcsec) 


960213 


1996Febl4 


50127.579 


L 


1.5"x6.0" 


3100-7200 


6 


2 


960214 


1996Febl4 


50128.432 


L 


1.5"x6.0" 


3100-7200 


6 


3 


960319 


1996Marl9 


50162.382 


L 


2.0"x6.0" 


3600-7400 


8 


3 


960321 


1996Mar21 


50164.395 


L 


2.0"x6.0" 


3600-7400 


8 


3 


960710 


1996Jull0 


50275.283 


L 


2.0"x6.0" 


3600-7000 


8 


1.4 


980120 


1998Jan21 


50834.632 


L 


2.0"x6.0" 


3800-7600 


8 


3 


980222 


1998Feb23 


50867.535 


L 


2.0"x6.0" 


3800-7600 


8 


2 


980504 


1998May04 


50938.293 


L 


2.0"x6.0" 


3700-7600 


8 


2.5 


980508 


1998May08 


50942.461 


L 


2.0"x6.0" 


3700-7700 


8 


2 


980725 


1998Jul26 


51020.703 


GH 


2.5"x6.0" 


3960-7231 


15 


2.2 


980726 


1998Jul27 


51021.712 


GH 


2.5"x6.0" 


3930-7219 


15 


2.5 


200126 


2000Jan27 


51570.959 


GH 


2.5"x6.0" 


4070-7350 


15 


2.5 


210513 


2001Mayl4 


52043.859 


GH 


2.5"x6.0" 


3980-7300 


15 


1.8 


210613 


2001Junl4 


52074.772 


GH 


2.5"x6.0" 


4022-7330 


15 


1.8 


210614 


2001Junl5 


52075.738 


GH 


2.5"x6.0" 


4010-7330 


15 


1.8 


21-07 


2001Jull6 


52107.320 


L 


2.0"x6.0" 


3630-8050 


8 


2 


220304 


2002Mar05 


52338.780 


GH 


2.5"x6.0" 


3976-7305 


15 


2 


220404 


2002Apr05 


52369.850 


GH 


2.5"x6.0" 


3976-7305 


15 


2 


220405 


2002Apr06 


52370.780 


GH 


2.5"x6.0" 


3976-7305 


15 


2 


220624 


2002Jun24 


52450.420 


L 


2.0"x6.0" 


3800-7600 


8 


2 


230126 


2003Jan27 


52666.96 


GH 


2.5"x6.0" 


3976-7305 


15 


2.5 


230127 


2003Jan28 


52667.930 


GH 


2.5"x6.0" 


3976-7305 


15 


2.5 


230325 


2003Mar26 


52724.900 


GH 


2.5"x6.0" 


3976-7305 


15 


3.5 


230522 


2003May23 


52782.750 


GH 


2.5"x6.0" 


3976-7305 


15 


3.5 
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Fig. 4. The temperature variation from 1996 till 2004. 

40%, as it was given in Popovic (2003). But, since we are 
investigating the changes in the BP during an interval the 
intrinsic reddening can be neglected as it should not vary 
too much in a relatively short period of around 8 years. 

4-2. Velocity Measurements 

These measurements were performed in order to present 
A as a function of Full Width Half Maximum (FWHM) 
and Full Width Zero Intensity (FWZI). The graphs A vs. 
FWHM and FWZI are useful because of: 

(i) the Balmer line ratios are velocity dependent in AGN 
(Stirpe 1990,1991) and this may be related to the physical 
conditions (temperature, density and kinematics) as well 
as to the radiative transfer effects; 

(ii) illustration of the physical conditions in the BLR 
plasma. As it was pointed out by Popovic (2003) if 
A < 0.3, the recombination "Case B" plus the intrinsic 



reddening may explain BP. The "Case B" recombination 
of Balmer lines can bring the log(F,i) vs Eu as a linear 
decreasing function (see Table 4.4 in Osterbrock 1989). 

In order to compare the temperature parameter A with 
the velocity of gas, very high data quality is required. The 
data quality is different for different spectra of NGC 5548. 
Consequently, for velocity measurements we used only the 
Ha and H/? lines. We first normalized the broad profiles 
of these lines and converted the wavelength axis to the 
velocity scale AT = (A — Ao)/Ao (the same as was given 
in Popovic 2003). We measured FWHM and FWZI of 
both lines, and we calculated the average values for each 
spectrum. 

5. Results and discussion 

5.1. BP of NGC 5548 

The results of our investigations are presented in 
Figures 2-4, and in Table 2. 

One can see in Figures 2 and 3 that BPs graphics indi- 
cate that the population of the excited levels from the 
Balmer series follows the Saha-Boltzmann distribution. 
We found that in all considered periods, the temperature 
parameter is A > 0.3. It means that the recombination 
"Case B" plus an intrinsic reddening cannot explain the 
line flux ratio of the Balmer lines. Also, we investigate 
the dependence of the temperature coefficient as a func- 
tion of FWHM and FWZI (Figure 3). As one can see 
from Figure 3, there is no significant correlation between 
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Table 2. The measured flux of the Balmer lines of NGC 5548, temperature parameter A. 



Spectra 




PHa/F 


H(3 






HP 




Phs/Fhp 




Phs/Fhp 


(erg- 


cm "^s 






A 




960213 


3 


405±0 


406 





396±0 


047 





174±0 


021 





071±0 


027 


(7 


998±0 


950)E- 


13 





456±0 


095 


960214 


3 


689±0 


625 





456±0 


074 





148±0 


025 





091±0 


016 


(7 


598±1 


220)E- 


13 





450±0 


085 


960319 


3 


647±0 


492 





341±0 


051 





133±0 


018 





077±0 


015 


(7 


281±0 


983)E- 


13 





515±0 


068 


960321 


3 


126±0 


384 





386±0 


048 





187±0 


029 





072±0 


016 


(7 


055±0 


864)E- 


13 





414±0 


100 


960710 


4 


203±0 


569 





465±0 


063 





204±0 


027 





084±0 


Oil 


(9 


891±1 


322)E- 


13 





462±0 


082 


980120 


3 


672±0 


446 





393±0 


048 





220±0 


038 





062±0 


008 


(8 


420±1 


023)E- 


13 





474±0 


138 


980222 


3 


624±0 


429 





376±0 


072 





210±0 


027 





082±0 


023 


(9 


694±1 


146)E- 


13 





426±0 


075 


980504 


3 


523±0 


398 





401±0 


048 





201±0 


024 





125±0 


042 


(9 


650±1 


090)E- 


13 





337±0 


013 


980508 


3 


571±0 


427 





416±0 


055 





212±0 


025 





117±0 


035 


(7 


998±0 


950)E- 


13 





344±0 


014 


980725 


3 


961±0 


471 





417±0 


049 





167±0 


020 





072±0 


009 


(9 


969±1 


179)E- 


13 





508±0 


088 


980726 


3 


727±0 


438 





338±0 


047 





153±0 


018 





039±0 


012 


(9 


235±1 


083)E- 


13 





633±0 


177 


200126 


4 


223±0 


528 





420±0 


052 





194±0 


027 




_ 




(7 


167±0 


886)E- 


13 





418±0 


022 


210513 


5 


718±1 


134 





314±0 


099 





109±0 


049 





146±0 


096 


(2 


900±0 


437)E- 


13 





586±0 


161 


210613 


4 


493±0 


622 





414±0 


056 





150±0 


077 




_ 




(3 


874±0 


524)E- 


13 





512±0 


073 


210614 


4 


439±0 


664 





330±0 
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widths and A. That is expected since macroscopic bulk 
motions arc mainly constrained by kinematics, probably 
an accretion disc as it was suggested by Shapovalova et al. 
(2004). On the other hand, it may indicate that changes 
in the self-absorption in the Balmer lines has no signifi- 
cant variation during the period. The averaged value ob- 
tained from all periods (full triangle in Figure 3) is in good 
agreement with the trend of the AGN sample given in 
Popovic (2003). Such a trend (increasing A with FWHM 
and FWZI) may indicate that physical conditions are de- 
pendent on kinematical characteristics of the BLR. This 
may be connected with distances from the massive black 
hole, e.g. the broader Balmer lines indicate higher macro- 
scopic bulk motions which can be expected in the regions 
closer to the massive black hole. 

Also, we could estimate the temperature using Eq. (3) 
& (4) from: 

T^^m. (12) 

In Figure 4 we plotted the temperature variation in the 
considered period: as one can see the averaged tempera- 
ture is wlOOOO K, and is changing by around 50% in the 
considered period. The maximum value of the temper- 
ature was reached in 1998 and minimum in 2002. The 
obtained temperatures agree well with the ones expected 
in the BLR. 

On the other hand, from the BP fits, we obtained the 
parameter B. As it can be seen from Eqs. (l)-(4) the 
parameter B depends on the number density of radiat- 
ing species and the partition function {B ~ \og{hcNo / Z)) . 



Both values are the same for the Balmer line series at one 
moment, but they are changing during the time, since we 
expect that the physical properties in the emitting BLR 
plasma of NGC 5548 are changing (seen as variability in 
the broad spectral lines). In principle the partition func- 
tion depends on the temperature and the level configura- 
tion (in our case Z is given for n ~ 2). It is hard to find 
some real information about variation in the number den- 
sity in the BLR plasma using the B parameter, but still 
one can have impression about magnitude of variations of 
these two quantities. Taking into account Eqs. (1) - (4) 
and also the fact that we normalized Balmer line fluxes 
to the H/3 flux, we can write log(A^o/^) ~ Be, -l-log(Fff^) 
[h and c are constants, and Bi is intercept corresponding 
to the energy of the lower level, since Z is given for n = 2, 
El — 10.2 eV). Using this relation we found the difference 
between maximal and minimal value of log(iVo/^) to be 
Alog(iVo/^) = log{No/Z),nax - log(A^o/^)mm ~ 0.9. This 
can be expected, e.g. if we roughly approximate that Z 
stay constant, then maximal variation of number density 
is 8 times, this is in an agreement with observed mag- 
nitude of variation in spectral lines of NGC 5548 (e.g. the 
flux of the H/3 was changing ~ 5 times.) 

5.2. Continuum vs. Temperature 

Here we start from the fact that the response of the 
broad lines to continuum variations in NGC 5548 suggests 
that BLR is very close to the continuum source (Ferland 
et al. 1992), or it may even indicate that a part of op- 
tical continuum is originating in the BLR as well. We 
assume that the continuum is originating in an accretion 
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Fig. 5. The continuum flux of NGC 5548 in minimum before 
(top solid line) and after (bottom solid line) subtraction of the 
host galaxy continuum (dashed line). 

disc, e.g. Kong et al. (2004) explained the variation in 
the optical spectral index with the variation in accretion 
parameters of the inner part of the disc. On the other 
hand, the Balmer lines may also be emitted from outer 
part of the accretion disc (Shapovalova et al. 2004), then 
the temperature variation measured from BP should be 
connected with the continuum variation. Therefore, we 
investigate the continuum intensity as a function of ob- 
tained temperatures. In an accretion disc, temperature 
depends on the radius (see e.g. Shakura & Sunyaev 1973) 
but if the continuum is emitted mostly from inner parts 
of the disc and Balmer lines mostly from outer parts, one 
can approximate that Tc correlates with the BLR temper- 
ature, i.e. Tc ^ Tblr. 

To find the correlation, we measured the flux in a win- 
dow from 4240 A to 4260 A (rest-frame wavelengths) and 
calculated an average value for 4250 A as well as in the 
spectral range from 5090 A to 5110 A (rest wavelengths) 
with an average value at 5100 A. 

In order to have only the AGN-component of the optical 
continuum, one should subtract the continuum of the host 
galaxy from the observed continuum of NGC 5548 nucleus. 
For this purpose, the contribution of the host galaxy at 
wavelengths 4250A and 5100A was obtained as following: 

a) according to the paper by Romanishin et al. (1995, 
Fig. 5) the host galaxy contribution in the AGN contin- 
uum of NGC 5548 at 5100 A in the aperture 5" x 7.5" is 
F(5100) = (2.9±0.54) • lO'^^erg • cm'^ • sec'^ • k-\ Since 
our spectra were taken with the smaller aperture 2.5" x 6", 
it is necessary to determine the correction for the aperture 
effect. We compared the continuum fluxes of NGC 5548 
at 5100A using spectra obtained in the same nights with 
our aperture 2.5" x 6" and with the aperture 5" x 7.5" of 
Peterson et al. (2002). The spectral observations dur- 
ing 13 different nights were used. The average correction 
for aperture effect is AF = F((5" x 7.5") - i^(2.5" x 6") = 
0.35 • lO^^^erg • cm~^ • sec~^ • A~^. The host galaxy con- 
tribution at 5100A with aperture 2.5" x 6" is F(5100)gal 
= (2.55 ± 0.54) • IQ-^^erg ■ cm~^ ■ sec^^ • A^^; 

b) For the determination of the host galaxy contribu- 



tion to the continuum flux of NGC 5548 at the wavelength 
4250A the spectrum of the normal E galaxy NGC 4339 
was used. The spectral observations of NGC 4339 were 
taken with the 2.1 m telescope with the same aperture 
2.5" X 6" and with the same spectral resolution (« 8 A) as 
the spectrum of NGC 5548 in the minimum activity state. 
After that we compared the spectra of NGC 4339 with the 
spectrum NGC 5548 in the minimum state (June 4, 2002). 
Such a comparison is justified because the central colors 
of the host galaxy NGC 5548 (B-VwO.9; V-Iwl.2), are 
similar to the colors of an elliptical galaxy (Romanishin 
et al. 1995). We varied the contribution from NGC 4339 
at 5100A from 50% to 100% of that in the observed spec- 
trum of the NGC 5548 and found, using a good substrac- 
tion of the absorption lines, that in the minimum activity 
of NGC 5548 the host galaxy contributes around 70% to 
the total continuum flux at 5100A (see Figure 5). From 
Shapovalova et al. (2004) the observed flux at 5100A 
in NGC 5548 for the minimum spectrum in aperture 
2.5" X 6" is F(5100)« 3.7 • lO^i^erg • cm'^ • sec"! • A"!; 
and then the host galaxy of NGC 5548 flux is F(5100)« 
(3.7 • 0.7) ■ 10-15 = 2.5 • lO-i^erg • cm'^ • sec"! • k-\ that 
coincides with the value obtained by us, using the data 
from Romanishin et al. (1995). 

We found that for NGC 4339 the fluxes ratio at 5100A 
and 4250A is F(5100)/F(4250)=1.45. Supposing that 
the continuum of the host galaxy NGC 5548 is similar 
to that in NGC 4339 we found that the flux of host 
galaxy NGC 5548 at A4250A is F(4250)=F(5100)/1.45 = 
(1.76 ± 0.54) • 10-i5erg-cm-2 -sec-i -A-i; Then we ob- 
tained the AGN-component continuum flux of the NGC 
5548 at A4250A and 5100A by substraction of the contin- 
uum flux of the host galaxy. The measured fluxes only of 
the AGN-component are given in Table 1. 

Note here that Peterson et al. (1995) showed that in 
the case of NGC 5548 the NLR has the same surface- 
brightness distribution as the PSF (i.e the NLR is point- 
like source ^2"), and that the point-source correction 
factor is always equal to one. Consequently, the broad- 
line/narrow-line flux ratio should be independent of the 
aperture and the seeing, at last for (slit)> 2", that is the 
case in our sample of observations. 

The total continuum flux (.^((A)) emitted at A can be 
represented as a sum of the BLR {Fblr), central source 
(Fc(A)) and stellar continuum (Fs(A)). As we noted 
above, the stellar continuum is subtracted and then: 

Ft{\T)^F,{\T,) + FBLR{X,T^w). (13) 

Assuming that the optical continuum is originating far 
enough from the centre of an accretion disc, one can use 
the next relation for the central source of continuum (see 
Peterson 2004, Eq. 4.6 and corresponding discussion) 

y8/3 

F,(A,T) ~ const X (14) 

If we neglect the contribution of the BLR continuum, 
the continuum flux ratio measured at Ai and A2 is 

Ft{\l,T) ,A2.i/3 . (^r.\ 
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There is no evidence that the continuum temperature 
is proportional to the BLR temperature, but one can as- 
sume that they arc in correlation i.e. that Tc ~ Tblr, 
therefore wc will present the continuum intensity vs. the 
BLR temperatures. 

In Figure 6, we present the AGN-component continuum 
flux variation as a function of the measured temperatures 
(at A4250A in Figure 6a and at 5100A in Figure 6b) and 
the flux ratio variation (Figure 6c). As one can see from 
Figure 6a, b, the AGN-component continuum tends to be 
a linear function of the temperature. We found a high 
level of correlation between the continuum and tempera- 
ture variability (r = 0.85). Also, we fitted the measured 
data with the continuum fimction as given by Eq. (14) 
and presented it with dashed line in Figure 6a, b, finding 
that Eq. (14) fits well the observed data. 

Furthermore, as one can see from Figure 6c the AGN- 
component continuum intensity ratio at 4250 and 5100 
A follows the expected function given by Eq. (15), i.e. 
remains constant as a function of the temperature. This 
also indicates the disc emission. But, here we should note 
that the temperature measured from the BLR fits very 
well the physical changes in the continuum. 

This can be expected because of the following reason: 
whether illuminated by a variable hard X-ray source at the 
centre, or by some other mechanisms, the optical/UV con- 
tinuum varies with the ionizing spectrum. Then, one can 
observe the continuum becoming "bluer" as it becomes 
brighter. This is consistent with an increase in the tem- 
perature of the thermal accretion disc as the central source 
(X-rays and the ionizing continuum emitting inner disc) 
brightens. The obtained high level of correlation between 
the AGN-component continuum and temperature, as well 
as of the ratio of the AGN-component continuum at dif- 
ferent wavelengths and temperatures, support a dominant 
emission by an accretion disc in the BLR. As it was noted 
by Uhich (2000) for the case of NGC 4151 the high cor- 
relation between X-ray emission (originated in an accre- 
tion disc) and temperature of BLR is expected. Here, we 
showed for the first time that this correlation exist in the 
optical part of the continuum. 

On the other hand, the correlation obtained from our 
measurements may also be explained by the following: the 
emitting gas in the BLR becomes hotter when the cen- 
tral source brightness and the flux of incident photons 
increases. At the same time the sizes of the Stromgren 
zones increase in depth, increasing the optical depths in 
the Balmcr lines. Their emitting efficiencies thus dimin- 
ish with increasing incident flux, but the effect on Ha is 
greater than that on H/3, which is greater than the effect 
on H7 (because Ha optical depth remains greater than 
H/3, etc., e.g. Ferland et al. 1979). In this case one can 
expect that the ratio of Balmer lines fluxes is a function 
of the continuum ratio measured in the blue and red part 
of the spectrum. To check it we presented the flux ra- 
tios of the continuum measured at 4250A and 5100 A as 
a function of Ha/ Hp, Ha/Hj and Ha/HS flux ratios 
in Figure 7. As one can see from Figure 7, a slight ten- 
dency exists for the line ratios F(Ha)/F(Hi) to be smaller 
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Temperature (K) 
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Fig. 6. The continuum flux measured at 4250 A and 5100 A 
as a function of the temperature measured with the BP (first 
two panels), and the ratio /(;7{4250)//(7(5100) as a function of 
the temperature (bottom). The measured values are denoted 
with full circles. The assumed linear function and the function 
given by Eq. (14) are presented with full and dashed line, 
respectively (first two panels). In the 3rd panel the full line 
represents the best fit with Eq. (15). 
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Fig. 7. The ratio of measured continuum flux as a function 
of the hne flux ratios Ha/Hfi (open circles), Ha/H^ (open 
triangles), and Ha/HS (open squares). In the bottom panel, 
the case for Ha/H6 flux ratio is presented separately. 

when the continuum is bluer, but it is a very tiny effect. 
It indicates that we have the real correlation between the 
temperature and AGN-componcnt continuum flux, that is 
presented in Figure 6. 

6. Conclusions 

By applying the BP method, as it was proposed by 
Popovic (2003, 2006ab), to the broad Balmer Hnes of NGC 
5548 observed from 1996 to 2004, we found that: 

(i) it seems that the coUisional processes play a signif- 
icant role and that the distribution of the excited level 
population (n > 2) follows the Saha-Boltzman equation, 
and this should be taken into account when modeling the 
BLR of NGC 5548; 

(ii) the BLR temperature was changed from 5000 K (in 
2002) to 15000 K (in 1998). The average temperature 
is 10000 K, which is a value expected for the BLR. The 
maximum of the reached temperature corresponds to the 
period where the Balmer lines were the most intensive; 

(iii) we found a correlation (r = 0.85) between the vari- 
ation of optical AGN-component continuum and temper- 
ature; this is the first time that this correlation is con- 
firmed, and it indicates the presence of an accretion disc 
in the BLR of NGC 5548 as it was earlier suggested by 
Shapovalova et al. (2004). 
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